Noise through Quantum Pumps 



(N 
O 

o 

(N 
X> 

IX, 



M. L. Polianski , M. G. Vavilov b < a , and P. W. Brouwer a 
a Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, NY 14853 
b Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455 

(February 1, 2008) 

We study the current noise through an unbiased quantum electron pump and its mesoscopic 
fluctuations for arbitrary temperatures and beyond the bilinear response. In the bilinear regime, we 
find the full distributions of the noise power and the current-to-noise ratio for a chaotic quantum 
dots with single-channel and many-channel ballistic point contacts. For a dot with many-channel 
point contacts we also calculate the ensemble-averaged noise at arbitrary temperature and pumping 
strength. In the limit of strong pumping, a new temperature scale appears that corresponds to 
the broadening of the electron distribution function in the dot as a result of the time-dependent 
perturbations. 
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PACS numbers: 73.23.-b, 72.10.Bg, 72.70.+m 
I. INTRODUCTION 

A periodic perturbation of a confined electron system 
may produce a direct current. The initial theoretical 
proposals p]|| and experimental realizations ^-^] for 
electron pumps were for a pump where the spectrum 
is gapped and the charge pumped in one cycle quan- 
tized and not subject to fluctuations. Recently, attention 
has shifted to pumps that are well connected to electron 
reservoirs, and, hence, do not have a gapped excitation 
spectrum [[|-^|. If the pump relies on a time-dependent 
perturbation that mainly affects the quantummechanical 
phases of the electrons, and not their classical trajecto- 
ries, it is referred to as a "quantum electron pump" . Such 
a quantum pump was fabricated by Switkes et al. [ fTo[ . 

The current that is pumped through a quantum elec- 
tron pump is subject to mesoscopic fluctuations and to 
quantum or thermal fluctuations (noise). Mesoscopic 
fluctuations of the current refer to the fact that the mag- 
nitude and direction of the time-averaged current vary 
from sample to sample. For a quantum pump built from 
a chaotic quantum dot, as is the case in the experiment of 
Rcf. pj], the mesoscopic current fluctuations were inves- 
tigated for various regimes of temperature, pumping am- 
plitude and dot conductances ]c^,p|,|lT|-p^| . On the other 
hand, noise — quantum and thermal fluctuations of the 
current — is a property of the current pumped through 
a particular realization of an electron pump. Noise in an 
electron pump is best described by the fluctuations of the 
charge pumped through the system in a certain number 
of pumping cycles. The statistics of such charge fluctu- 
ations was studied in Refs. |l5|-|2^] for temperatures and 
pumping frequencies much smaller than the inverse dwell 
time (escape rate) of electrons in the quantum dot. 

In this paper, we consider the mesoscopic fluctuations 
of the noise. We do not impose any restrictions on the rel- 
ative magnitudes of temperature T, pumping frequency 
u>, escape rate 7 and pumping amplitude. This is im- 
portant, as in an experimental realization of a quantum 



pump, T and 7 are usually comparable, while both are 
much larger than u>. Previous works by Andreev and 
Kamenev [jl5| and Levitov |l6| addressed the full count- 
ing statistics, but at temperatures kT -C hid only. The 
mean square charge fluctuations for hto, kT <C 7 (but in- 
cluding the case tko kT) were considered very recently 
by Moskalets and Biittiker 

Denoting the quantummechanical average with a bar 
TT7 , the quantum and thermal fluctuations of the pumped 
charge are described by 
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Here To is the observation time and Q is the total charge 
pumped through the dot in the time tq. The noise in an 
electron pump can be divided onto a Nyquist-Johnson 
component S N and pumping component S p . The for- 
mer is the thermal equilibrium noise due to the thermal 
fluctuations of electrons in the leads, and depends on 
the electron temperature T in the leads and the time- 
averaged conductance G of the quantum dot through the 
fluctuation-dissipation theorem, 
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The pumping contribution S p , in contrast, is a true non- 
equilibrium noise due to the perturbation of electrons 
inside the dot. 



Ill 



this con- 



As we shall discuss in Sec. 
tribution to the noise can be seen as arising due to the 
heating of electrons inside the dot as a result of the time- 
dependent perturbation. 

In the adiabatic regime hto <C 7 one needs to vary at 
least two system parameters periodically in order to gen- 
erate a direct current. Current noise, however, is already 
generated if only one parameter is varied. The problem 
of current noise (and full counting statistics) for a sin- 
gle time-dependent scatterer with frequency Hlo <C 7 was 
addressed by Levitov and coworkers |n||2(J. Our results 
can be used to compute the mesoscopic fluctuations of 
the current in that case. We find that the main effect 
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of the second time-dependent parameter in a true elec- 
tron pump is to reduce the mesoscopic fluctuations of the 
noise. 

The paper is organized as follows: A formal expression 
for the noise in terms of the time-dependent scattering 
matrix of the dot S is derived in Sec. |l| Section |ni] con- 



siders the mesoscopic fluctuations of the noise through a 
chaotic quantum dot with two (or more) time-dependent 
parameters. For an adiabatic quantum pump at tem- 
perature kT <C 7 we consider the full distribution of 
the mesoscopic fluctuations of the noise in Sec. Ill A. 



We focus on the cases of a quantum dot with single- 
channel and many-channel point contacts. For a dot with 
many-channel point contacts, the sample-to-sample fluc- 
tuations of the noise are much smaller than the average. 



In Sec. [II B we then present the ensemble-averaged noise 



for a dot with many-channel point contacts at arbitrary 
temperature and pumping strength. We conclude in Sec. 



IV . Finally, in Appendix we give some details of the 
calculations of the averages over the ensemble of chaotic 
quantum dots. 



bations, so that the system is described by a scattering 
matrix that depends on energy, but not on time. On the 
other hand, as long as pumping frequencies and tempera- 
tures much smaller than the escape rate 7 are considered, 
the system with a time-dependent perturbation can be 
described using a scattering matrix that depends on time, 
but not on energy |I^ , |l6| , |l9| -p2"| . When kT, Tvuj and 7 are 
all comparable, one needs to use a scattering matrix S 
that depends on both time and energy, or, equivalently, 
that depends on two times or two energies. Here we 
shall make use of a formulation with a scattering matrix 
S(t, t') that depends on two times. This formulation was 
used to calculate the the time-averaged conductance and 
pumped current in the Ref. |l3f| . The formalism is equiv- 
alent to the two-energy scattering formalism developed 
by Biittiker and coworkers for time-dependent transport 
through mesoscopic structures |2^-|27|| . 

The scattering matrix S(t, t') relates the annihilation 
[creation] operators a Q (t) [a T (f)] and b Q (t) [bt(i)] of 
incoming states and outgoing states in channel a — 
1, . . . , 2N of the leads (the index a includes the spin de- 
gree of freedom) , 



II. GENERAL FORMALISM 

We consider fluctuations of charge transmitted the 
sample during the observation time To. The sample is 
connected to electron reservoirs through two point con- 
tacts with N L and Nr open channels, respectively, see 
Fig. |l} It is subject to a periodic perturbation at fre- 
quency uj/2ir; The perturbation is described by speci- 
fying the time-dependence of parameters Xi(t), Xzft), 
. . . , X n (t) characterizing Hamiltonian of the sample. The 
electrons in the two reservoirs are held at the same chem- 
ical potential fi and temperature T at all times during the 
pumping cycle. 
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FIG. 1. Schematic drawing of the quantum dot and the 
leads. There are Nl {Nr) propagating modes in the left 
(right) contact. The shape of the dot is controlled by the 
voltages of two shape-defining gates. The vectors bll,r and 
^>l.r of annihilation operators for the incoming and outgo- 
ing states, respectively, in the left (L) and right (R) leads are 
related by the scattering matrix 5. 

We calculate the noise of the quantum pump using 
the scattering formalism of Biittiker p3[ . In Buttiker's 
original application there are no time-dependent pertur- 
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Here the indices a and (3 label the propagating chan- 
nels in the point contact contacting the dot to the left 
and right reservoirs for a, (3 = 1,...,2Nl and a,0 = 
2N L + 1, . . . , 2N, respectively. The 2N x 2N matrices 5+ 
and S are related as 



(5t(f,t)) =S* 0a (t,t'). 



(4) 



Causality requires that S(t,t') = if t < t'. We restrict 
our attention to the case where spin rotation invariancc 
is preserved, and, hence, S is proportional to the 2x2 
unit matrix in spin grading. 

The expression for the current I L in the left lead is 



2N L 



lL(t) = eY] [at (t)a„(t) - bt (t)b„(t)] 



a=l 



(5) 



A similar expression holds for the current 1^ in the right 
lead. Since we are only interested in the charge pumped 
through the dot in a long time interval and the charge 
on the dot is conserved after each cycle, we replace the 
expression for the current operator by a suitable combi- 
nation of Ii and Ir , 

I(t) = ^I L (t) - ^I R (t) 

2N 

= e (4(*)A«0B0(t) " *>t (tJA^Mt)) . (6a) 

a.fi=l 
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Here A is a diagonal matrix with elements 
' N R /N a = l,...,2N L , 



-N L /N a = 2N L + l, 



,2N. 



(6b) 



In terms of the current operator I, the time-averaged 
pumped current I reads 



1 f T ° - 
I = — dtl(t), 
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where the observation time To is an integer number of 
pumping cycles. The noise S is defined as [cf. Eq. (Q)] 

(8) 
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S = — I dtdt' I(t)I(t') - I(t)Ifi') 
to Jo v 

In the leads, the electron distribution function is given 
by the Fourier transform f(t) of the Fermi function, 



al(t')a a {t)=8 a pf{t'-t), 
SL f} (t')f&(t) = 6 afi f(t-t), 



(9) 
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where we defined 



f(t) = 6(t)-f(t). 



(10) 



Using Eq. (H) to eliminate the operators b(t) and b^(t) 
from Eq. ( pa] ) and Eq. ^ to compute the quantumme- 
chanical expectation values, we find the average current 

I = —e [ dt [ dtidt 2 f(ti - t 2 ) 
to Jo J 

xTr [S t (t 1 ,t)AS(t,t 2 )-S(t-t 1 )A5{t-t 2 )] 1 (11) 
and the noise 

S = S N + S P , (12) 

where the Nyquist- Johnson and pumping contributions 
S N and S p read 



S N = —J dtdt' I dhdt 2 f(ti -*')/(*'- *a) Tr [S(t - t 1 )A 2 5(t - t 2 ) - S f (<i, t)AS(t, t 2 )A] , 

P 2 /-to /• 

S p = — dtdt' / dhdtidt^dt'zfit! - t' 2 )f{t[ - t 2 ) 

TQ Jq J 

x Tr [«S%,t)A<S(t, ^2)5+ (*;,*') A5(t', t' 2 ) - <5(t - - t[)A 2 S(t - t 2 )5{t' - t' 2 )] . 



(13) 



(14) 



(These equations can also be derived using the Keldysh 
formalism, see Ref. p8[ .) 

For adiabatic pumping, Hu) <C 7, Equation rtllj) is 
equivalent to the time-averaged current of Refs. |S|] and 
H . The Nyquist- Johnson contribution to the noise is re- 
lated to the time-averaged conductance G of the system 
at temperature T, see Refs. and Eq. (0) above. 



A. Bilinear adiabatic pumping 

Of particular interest is the case when the perturba- 
tion is slow compared to the (elastic) escape rate 7 of the 
electrons from the sample into the reservoirs. This is the 
regime of the adiabatic quantum pump of Ref. fj(i|] . In 
this approximation it is advantageous to use an analog 



of the Wigner transform for the matrix S(t, t'): 

/+00 
dt'e-^-^Sfrt'). (15) 
-00 

Up to corrections of order tkj/j, the matrix S(e, t) is 
equal to the "instantaneous" scattering matrix Sx(e), 
which is obtained by "freezing" all parameters Xj to their 
values at time t pj]J3^ |. 

If, in addition to being adiabatic, the parameters Xj, 
j = 1, . . . ,n, undergo only small excursions from their 
average value, which we set to zero, we may further ex- 
pand in Xj . To find the noise it is sufficient to expand S 
up to the second order in X. Arranging the parameters 
Xj in an n-component vector X = (A/, . . . , X n ) T , we 
thus find 



S p = — dtdt' 
4r Jo 
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x [(X T (t) - X T (t , ))JC(e,e')(X(t') - *(<))] , 
where /(e) = 1 — /(e) and the n x n matrix K. reads 

fCij{e,e') = Tr [A 2 K t {e)Kj{e) + A 2 ^(e')^ 2 (e') - 2kn i (e)Afc j {e')] , Kj{e) 



£ [ + °° ^fcos[(e -s')(t-t')} /(.-)./(.-') ^. /,.-)/( r'l 
27r Zir 



■ dSx(e) ot 



dXi 



(16a) 



(16b) 
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In the bilinear regime and at zero temperature, the 
time-averaged current I and the noise S p are suffi- 
cient to define the total counting statistics of the pump 
JUsf : To lowest order in the excursions of the param- 
eters Xj, the pumping cycles are statistically indepen- 
dent. A pumping cycle is characterized by quantum- 
mechanical probabilities Pr and Pj, that an electron is 
pumped from left to right or from right to left, respec- 
tively. Both Pa and Pl are small, of order X 2 , and one 
has I = (eu/2ir)(P L - P R ), S p = {e 2 uo/2n)(P L + P R ). 

For a typical quantum dot with capacitance C and 
mean level spacing A, the charging energy e 2 /2C 3> A. 
The derivatives in Eq. ( 16b ) should be taken at a constant 
value of the chemical potential /i, which, in the Hartree 
approximation, is equal to sum of the electron's kinetic 
energy and the electrostatic potential. In the absence of 
electron-electron interactions all derivatives are taken at 
constant value of the kinetic energy. As we prefer to take 
derivatives at constant kinetic energy (i.e., at constant e) 
in both cases, we substitute the parametric derivatives as 

MM. 



d 

dX 



d 

dX 



TTr 



de 



Equation ( )l7| ) differs from a similar expression in Ref. 
H by factor of 1/2 in front of the traces because of the 
double size of matrix S as a result of the inclusion of spin. 



and independently distributed. We choose the scale for 
the parameters Xi such that the mean square deriva- 
tive {(de^/dXiXds^/dXj)) = ^A 2 /^ 2 , where A is the 
mean level spacing and is an energy level in the closed 
dot. 

The transmitted charge Q is measured during a time 
To, which we will assume to be a large number of pump- 
ing cycles. This requirement of large observations is dis- 
cussed in detail in Ref. [^0) . For short observation times 
boundary effects related to switching processes in the sys- 
tem need to be taken into account. 



A. Weak adiabatic low temperature pumping 

If not only the frequency huj is much smaller than the 
escape rate 7, but also kT <C 7, the scattering matrix 
S(e,t) in Eq. ( [l6| ) can be taken at the Fermi level e = 0, 
and we find the simple result 



S p = — dtdt'fit - t')f(t' - t) 
to Ja 

x 1 [(X(t) T - X{t') T )lC{X{t') - X(t))] , (18a) 



Idj = Tr (A 2 K i K : j + A^TljUi ~ 2A1Z l A1l 3 ) 



(18b) 



We now consider the case of two time-dependent pa- 
rameters, X(t) = (Xi(t),X 2 {t)) T : 



III. APPLICATION TO CHAOTIC QUANTUM 
DOTS 

We now consider the mesoscopic fluctuations of the 
pumping noise S p for the case of a chaotic quantum 
dot. The quantum dot is characterized by a mean level 
spacing A, escape rate 7 = NA/2n, Thouless energy 
-Erh = fy/ferg 3* A, and capacitance C, with charging 
energy e 2 /2C ^> A. It is coupled to two electron reser- 
voirs via ballistic point contacts with Nl and N R chan- 
nels each, see Fig. [l| 

The dot is driven by periodically varying parameter(s) 
Xi, i = 1, . . . ,n with frequency u>, cf. Fig. 0. In the 
experiment of Ref. fTo[| , the parameters Xi correspond 
to the voltages on external gates that control the dot 
shape. The precise relation between the parameters X 
used in the theory and in the experiment is not known 
a priori, but can be established using independent mea- 
surements of, e.g., the derivative of the conductance |3^] 
or the rate of change of the position of Coulomb blockade 
peaks when the point contacts between the dot and the 
reservoirs are pinched off p3| . Following Refs. || , fl3[ , 
we will assume that the different parameters Xi corre- 
spond to different perturbations of which the matrix ele- 
ments between states (of the closed dot) within a Thou- 
less energy from the Fermi level e = are Gaussian 



Xi{t) = Xx cos(wt), X 2 (t) = X 2 cos(ujt + 4>), (19) 



in more detail. For this case, Eq. ( |l8a| ) is factorized into 
a factor F(T,u,to) that depends on the relevant time 
and energy scales, and a sample specific contribution, 



S p = e 2 F{T,iu,T ) {KuXf + K 22 X. 
+ 2cos0 IC 12 X 1 X 2 ) . 



(20) 



To find the integral F(T, uj,tq) for the physically relevant 
limit of long observation times tq ^ 1/T, we note that 
in this limit the distribution function f(t) is given by 
Fourier transform of the equilibrium Fermi distribution 
function, 



/(*) 



de 



iet/h 



ikT 



2irh e e / kT + 1 2% smh{nkTt/h) ' 



Substitution into Eq. (18a) then yields 

,_, uj ( . Tllu 2kT 

F = — coth 

2tt V 2fcT huj 



(21) 



(22) 



In the opposite limit tq <C Ti/kT, F{T, u>, tq) can be found 
using the procedure proposed by Ivanov and Levitov p0[ . 
In a finite observation time To, energy levels are resolved 
up to H/tq, so that we must set 
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/(*) 



1 



2r sin [ir(t + iO)/r ] 



(23) 



Then, following Ref. (2Cj], we change the time variable 
t e [0, To] to z = e 27Tlt rn>^ anc j perform contour integra- 
tion in the circle \z\ = 1 in complex plane. The result co- 
incides with the low-temperature limit of Eq. ( p2|) above. 

The physical meaning of the function F(T,uj,to) with 
To ^> 1/T becomes clear once it is written in energy rep- 
resentation, 



F 



ds 
1 



(24) 



Equation ( p4| ) measures the change in the number of equi- 
librium electron-hole pairs due to absorption and emis- 
sion of the pumping field quantum hto. At temperatures 
kT <C huj the Fermi distribution is sharp, and F — w/2-7r. 
At high kT ^S> fiui, the Fermi distribution is smooth on 
the scale ~ w, so that F is small, F ~ Twj 2 /kT. 

Fluctuations of the noise are described by the second 
factor in Eq. (po|), which varies from sample to sample. 
To find the mesoscopic fluctuations of the noise, we use 
the joint distribution of the matrices Hi (i = 1, 2) derived 
in Ref. |34|] . The resulting distribution of this mesoscopic 
contribution depends on two parameters, 



Xl), C c =-X 1 X 2 sin 



(25) 



Equality C c — C\ is achieved if X\ — X 2 and cf> — t/2, 
corresponding to the circular contour in the (Xi,X2) 
plane. Below we consider the noise distribution for a 
quantum dot with single-channel point contacts (Nl — 
Nr = 1), and for a dot with many-channel point contacts 
(Nl , Nr > 1) separately. 



1. Two channel geometry, N — 2 

For small N, the full distribution of the noise can be 
obtained using the method of Ref. |55| to numerically 
generate the matrices IZi and IZ2 according to the ap- 
propriate distribution, see Appendix for details. Fig- 
ure ^| shows the distribution of the noise power S p and 
the current-to-noise ratio I/S p for a quantum dot with 
single-channel leads (N — 2) for the case C\ = C c . The 
distributions are shown with and without time-reversal 
symmetry (TRS). For reference we have also included the 
distributions for the case e 2 /C < A of weak electron- 
electron interactions inside the dot. The case C c < C\ 
(i.e., the dependence on the phase difference cj>) is illus- 
trated in Fig. ^|. The current-to-noise ratio shown in 
Figs. |^ and |^ accounts for the pumping noise S p only. 
The Nyquist-Johnson noise S N presents a different noise 



source, which will dominate over the pumping noise if 
kT > Tiuy/Ci, 
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FIG. 2. Main panel: Distribution of pumping noise S p for 
a chaotic quantum dot with single-channel point contacts and 
two time-dependent parameters X\ and X2 given by Eq. ( p"s| ) 
with Ci = C c - [The parameters Ci and C c are defined in Eq. 
(p5[).] The noise is measured in units of e 2 CiF. The plots 
are with and without time-reversal symmetry (TRS) and for 
e 2 /C <C A (weak electron-electron interactions inside the dot) 
or e 2 /C 3> A. Inset: Distribution of the current-to-noise ratio 
I/S p , measured in units of w / (2nFe). There is no divergence 
at I/S p -> 0. 



We note that the distributions of S p and I / S p are 
highly non-Gaussian. In particular, the mean (S p ) of 
the noise distribution is dominated by the algebraic tail 
for large S p , and is not representative of the distribution 
itself. [For example, in the absence of time-reversal sym- 
metry and for e 2 /C < A, the mean (S p ) = 8/3(e 2 QF), 
while the most probable value is for S p « 0.5(e 2 C/F).] 
We also note that, while the phase difference (j> affects 
the typical size of the time-averaged pumped current 
/ oc ewC c but not the form of the distribution, chang- 
ing <fi has a small effect on the average noise (S p ), but 
changes the shape of the noise distribution significantly, 
see Fig. |[ In particular, the probability to find small S p 
is significantly higher for (j) close to zero than for ip ~ n/2. 
The reason for this difference is that the case <fi ~ tt/2 cor- 
responds to noise generated by two independent sources, 
while 4> close to zero corresponds to only a single noise 
source. Furthermore, the current-noise ratio I/S p has 
a maximum at I/S p — e, as was predicted by Levitov 
fujfl . For a quantum dot with single-channel point con- 
tacts, there is a finite probability density to achieve this 
optimum current-to-noise ratio, as is seen in the inset 
of Fig. ^|. For point contacts with more than one chan- 
nel, the probability density to attain the maximum value 
I / S p = e vanishes [16|, see, e.g. , Fig. |^. 
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FIG. 3. Distribution of the pumping noise S P (main 
panel) and the current-to-noise ratio I / S p (inset) for a quan- 
tum dot with two single channel point contacts and broken 
time-reversal symmetry, for various values of C c /Ci. The val- 
ues of Cc/Ci shown are C c /Ci — 1, sin(7r/4), sin(7r/12), and 0. 
If the variations of the parameters X\ and X2 have equal am- 
plitudes, this corresponds to phase difference cj> = 7r/2, 7r/4, 
7r/12, and 0, respectively. The noise S F is measured in units 
of e 2 CiF and the current I is measured in units of eCiu)/2n. 
For C c = 0, P(I/S P ) is a delta function at I/S p = (not 
shown) . 



2. Multichannel limit, N S> 1 

For large N, the ensemble average and variance of the 
noise can be expressed in terms of an integral over the 
unitary group and over the eigenvalues of the Wigner- 
Smith time-delay matrix K = -ifi{dS / de)S^ §§. These 
integrals can be calculated using the method of Refs. 
p6| , |37t , together with asymptotic expressions for the den- 
sity and two-point correlations of the eigenvalues of the 
Wigner-Smith time-delay matrix, see Appendix for de- 
tails. 

The result is 



(S p ) = 2e'Fg ( 1 + j^S 0A ) Q 



(26a) 
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1 + (1 + 2^,1) 
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(26b) 



where we abbreviated g — NlNr/N . The factor F de- 
pends on lu and T and was defined in Eq. (|22|). Note 
that the fluctuations of S p are a factor ~ 1/N smaller 
than the average. Higher cumulants of the noise are even 
smaller, so that we conclude that, for N ^> 1, the noise 
distribution becomes sharply peaked at the average (S p ) 



and that the remaining sample-to-sample fluctuations are 
Gaussian. 
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FIG. 4. The same as Fig. |^ for a quantum dot with 
Nl = Nr = 5. Curves with weak and strong electron-electron 
interaction are indistinguishable. 

In Fig. H we show the results of a numerical calcu- 
lation of the distribution noise S p and the current-to- 
noise ratio I/S p for a dot with N L = N R = 5, N = 10. 
This value of N can be seen as intermediate between 
the small- N regime, where the distributions are strongly 
non-Gaussian and the large- N regime where the distri- 
butions are Gaussian. We note that for N — 10 the noise 
distribution still has pronounced tails for large S p . 



B. Noise at arbitrary pumping for multichannel dots 

In this subsection we consider a general time depen- 
dence and amplitude of the parameters Xj. We limit our- 
selves to the calculation of the ensemble averaged noise, 
since the mesoscopic fluctuations of the noise are smaller 
than the average by a factor 1/N if N 1, see Eq. (|2rj). 

For a calculation of the ensemble average noise (S p ), 
we need to know the correlation functions of scatter- 
ing matrix elements S a p(t,t') for an ensemble of chaotic 
quantum dots. The indices a and refer to the "or- 
bital" channels as well as to the spin of the electrons. In 
order to discriminate between the orbital and spin de- 
grees of freedom, we set a = (i, a), where a = ±1 refers 
to spin and i = 1,...,N denotes the orbital channels, 
i = 1 , . . . , Nl for channels in the left point contact and 
i = Nl + 1, . . . , N for channels in the right point contact. 
In this notation, the scattering matrix S a /3 — S°j5acr' is 
proportional to the 2x2 unit matrix in spin space. The 
correlator of the orbital part reads 



(SM)S k t(T,r')) = 6(t-t'-T + r') 



t) 



, (27) 
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with T> and C given by 



T>(tx,t 2 ,T) =e(ti-t 2 )exp{- 



C(t u t 2 ,T) =6(t 1 -i 2 )exp. 



1 A J [7V d + 2 (X T (£ - r/2) - A- T (C + r/2))(*(£ - r/2) - *(£ + r/2))] | , (28) 



ti — r/2 



t 2 — r/2 



Ad?7 
2irh 



[N c + 2 (A- T (?7) - X T (h + t 2 - r]))(X(v) - X(h +t 2 - 7]))] 



(29) 



Here Q(z) 
ated 



1 if z > and otherwise and we abbrevi- 



N,, 



N, N c = N + 2x 2 , (30) 

where x oc ( < &/$o)(' r crgA)~ 1 / 2 is a dimensionless parame- 
ter describing the magnetic flux penetrating the quantum 
dot and $o the flux quantum p9| ]. The unitary ensemble, 
when time-reversal symmetry is fully broken, corresponds 
to the limit x — > oo. 

In the literature, two equivalent approaches have been 
taken to calculate correlators such as Eq. (U7\ ) above, 
the Hamiltonian and scattering approaches J39J . In the 
Hamiltonian approach the fundamental object is the ran- 
dom Hamiltonian of the closed chaotic quantum dot and 



the Green functions related to it. Once the scattering 
matrix S is expressed in terms of Green functions, the 
correlator ( |27j ) can be analyzed by standard diagram- 
matic techniques. The two terms C and T> then appear 
as the cooperon and diffuson contributions. The funda- 
mental object of the scattering approach is a statistical 
model for the scattering matrix S of an ensemble of dots. 
Equivalence of both methods, including the parametric 
and energy dependence of S is shown in Ref. j34|. A 
derivation of Eq. (27) using the scattering approach is 
given in Ref. 40 1. We refer to Ref. |^8| for a discussion 
based on the Hamiltonian approach. 

Knowing the correlator (|27j), we can find the ensemble- 
averaged noise (S p ) 7 



2e 2 g 

TO 



dt dt'f(t - t')f(t' - t) 



t + t' t + t' h( 

2 '~ 2 ~ 



,t'-t) d( 



(31) 



where g = NlNr/N is the dot conductance, 7 = NA/2tt 
the escape rate, and we assumed tq ^> l/u),h/j. The 
Cooperon contribution C of Eq. ( p7| ) does not contribute 
to (S p ) to leading order in 1/iV; its contribution is a 
factor 1/N smaller, see, e.g., Eq. (26a). 

Equation (|3l]) gives the ensemble-averaged noise for 
arbitrary to, T, and 7, and for arbitrary excursions of the 
parameters Xj . We now investigate Eq. (|l|) for the case 
that there are two time-dependent parameters X± and 
X 2 with time dependence given by Eq. (p^[). In order 
to distinguish regimes of "weak" (bilinear) and "strong" 
pumping, we introduce the energy scale []13|,|30[ 



kT* = huj. 



(32) 



The meaning of the energy scale kT* becomes clear when 
we view the pumping process as "diffusion in energy 
space" : carriers absorb or emit energy quanta of size hu) 
at a rate X 2 A/h [Q. Weak (bilinear) pumping corre- 
sponds to the case when the probability to absorb or emit 



one or more quanta is small, i.e., to the regime X 2 A <C 7 
or kT* -C hu). For kT* ^> hui many quanta are absorbed 
or emitted, so that the carriers in the dot shift their en- 
ergies by an amount ~ huj \/X 2 A/7 ~ T* . If T* exceeds 
the temperature T of the electrons in the leads, the time- 
dependent potentials in the dot lead to significant "heat- 
ing" of the electrons inside the dot and T* can be viewed 
as an effective electron temperature inside the dot. The 
latter regime is referred to as "strong" pumping. 

For weak pumping, kT* <C hu>, we can expand T> to 
first order in X 2 and we recover the result Eq. (26a), now 



without a restriction on the temperature T. For strong 
pumping, kT* 3> max{?ia;, kT}, a simple expression for 
the noise power can be obtained if pumping is adiabatic, 
Tiuj -C 7. In that case, we note that T> in Eq. ( |3l| ) 
contains the fast decay ~ exp(— £) and a slowly varying 
contribution from the time-dependence of the parameters 
Xj. Since the Xj vary slowly on the time scale %/'), the 
integration over £ can be done, and the result is 



TO JO 



dt'dtf(t' - t)f(t - t') 



N 



N + 2(X T (t) - X T (t'))(X(t) - X(t')) 



- 1 



(33) 



In the limit of low temperatures, kT <ti hu (and, as before, assuming long observation times h/ro <C kT,hco), Eq. 
© yields 



P . _ e 2 guj [1 + 6Q(1 + Q(l - S 2 ))} E(fc) - [1 + 2C/(1 - S)} K(fc) 



tt[1 + 2Cj(l - S)]y/1 + 2Q(1 + S) 



(34) 
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where E(fc) and K(fc) are full elliptic integrals of the 2nd and 3rd kind, respectively, and we abbreviated 



CcY 2 4Ci s 

CJ ' l + 2Cl(l + S)" 



(35) 



In the special case that the two time-dependent parameters X\(t) and -X^t) have equal amplitudes, S — | cos</>|. The 
dependence of the averaged noise in Eq. ( |34| ) on the ratio C c /C; (i.e., on the phase difference <j>) is weak. For the case 
C c = (0 = 0) wc find the asymptotes 



{S p )=e 2 g-C l= e 2 g 

7T 



(fcT*) 2 
■nh UJ 



if fcT* < Tiw 



7T^ 



/ C P\ 2 3w /C/ o 
(b ) = e q—^\ — = e o 



fcT* if fcT* > Tkt) 



(36a) 
(36b) 



For the case C c — 0, but at arbitrary temperatures fcT, Eq. (|33|) yields 



(S P ) = e 2 gco 



2fcT 



2(kT*) 2 hLO sin 2 C 



sinh^(27rfcTC/?iw) 



^/(M 2 + 4(fcT*) 2 sin 2 C [(M 2 + 4(fcT*) 2 sin 2 C] 



2 Al 3 / 2 



(37) 



The low temperature asymptotics can be easily obtained 
from this result, and reproduce the Eq. ( |36| ) for fcT* <C 
Tia; or fcT* ^> 7kj. For intermediate values of kT* jhu) and 
kT/huj, Eq. @ is plotted in Fig. § As long as both Tilj 
and fcT are much smaller than fcT*, the integral ( p^ ) 
is dominated by £ <§; TiujjkT* <C 1, so that the strong 



pumping asymptote of Eq. (36b) is reached, irrespective 
of cj or T. 




FIG. 5. Noise (S^) in units of e 2 gu>, as a function of 
the dimensionless pumping strength kT* /hui and for various 
choices kT/Tioo. 



An analytical expression for <p ^ can, in principle, be 
obtained from Eq. ([53|) as well. The qualitative behav- 
ior as a function of pumping strength and temperature is 
similar to that shown in Fig. || for the case cf> = 0. The 
limit of strong pumping should converge to the limit Eq. 
(531) of T = 0, similarly to the case = studied above. 



IV. CONCLUSION 



In this paper we calculated the current noise generated 
in a quantum pump using a scattering matrix formal- 
ism with a scattering matrix S(t, t') that depends on two 
times. With this formalism, we could consider arbitrary 
pumping frequency u>, temperature T, escape rate 7, and 
pumping strength X. We then calculated the average 
and variance of the noise for an ensemble of quantum 
pumps consisting of a chaotic quantum dot. 

One issue that has received considerable attention re- 
cently is the question whether one can build a noiseless 
quantum pump ^T| . While our results for the bilinear 
pumping regime show that there is a finite (mesoscopic) 
probability for zero noise, it is not possible to have no 
noise and a finite pumped current at the same time in 
the bilinear regime fll6|]42] ], cf. Figs. || and |j| This is 
different beyond the bilinear regime, where a quantized, 
and, hence, noiseless pump has been proposed using a 
pumping contour that encircles a resonance in an almost 
closed quantum dot pjjfl . 

We note that both the equilibrium Nyquist-Johnson 
noise S N and the pumping noise S p depend on the avail- 
able energy window. For pumping noise, that energy 
window is the heating temperature T*, see Eq. (|32|); for 
Nyquist-Johnson noise it is the temperature T of the 
electron reservoirs. The results (^) and (|34j) (^7|) allow 
us to compare the Nyquist-Johnson and pumping con- 
tributions to the averaged noise. In the experimentally 
relevant case that Tioj <C fcT, both noise contributions 
are proportional to the (dimensionless) dot conductance 
g, but the Nyquist-Johnson noise scales as T, while the 
pumping noise scales as T* if T* > T and as T* 2 /T if 
T* -C T. The Nyquist-Johnson and pumping contribu- 
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tions to the noise are comparable at T ~ T* . An experi- 
ment cannot separate the two contributions to the noise, 
since it measures the total noise power. The pumping 
contribution to the noise is dominant as long asT-CT*. 
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APPENDIX: INTEGRATION OVER THE 
MATRICES TZ. 

For the integration over the matrices 1Z\ and 7^.2 we 
make use of the fact that they can be parameterized as 



Kj = -i—Uf^HjT^U" 1 <g> l s 
zttH 



1,2, 



(Al) 



where U is an N x N unitary matrix, f is a diagonal NxN 
matrix containing the eigenvalues r m , m — 1, . . . , N, of 
the Wigner-Smith time-delay matrix on the diagonal [B8| , 
Hj is an N x TV" hermitian (real symmetric) matrix is 
time-reversal symmetry is broken (present), and I2 is 
the 2x2 unit matrix in spin grading. 

For a chaotic quantum dot, the distributions of the her- 
mitian matrices H\ and Hi, the unitary matrix U, and 
the diagonal matrix f are all independent. The matrices 



Hj, j = 1,2, have a Gaussian distribution, 



P(H) cx exp(-/3TrH 2 /8), 



(A2) 



where (3 = 1 if time-reversal symmetry is present and 
P = 2 if time-reversal symmetry is broken by a magnetic 
field. The matrix U is uniformly distributed in the uni- 
tary group, and the eigenvalues r m , m = 1, . . . , N of the 
time-delay matrix have distribution 



N 



N 



Pa l 



nhC 



eh in 

m— 1 / m<n 



1 m 1 n I 



N 



- .--3PN/2+P-2 -Pirn/A-r 

1 rn. 1 > m c 



(A3) 



m— 1 



Knowing these distributions, finding the distribution 
for small N becomes a matter of mere quadrature. We 
have obtained the plots of the distributions of S p and 



I/S p by numerically generating 10 7 -10 8 matrices TZj dis- 
tributed according to the above distribution. We refer to 
Ref. |35| for the details of implementation of this proce- 
dure. Moments of the noise and the current can be found 
by performing the Gaussian integrations over H and the 
integrations over the unitary group with the help of the 
technique of Ref. |36|. For small N, the remaining inte- 
gration over the r m can be done explicitly. For large N, 
it is sufficient to know the density and two-point corre- 
lator of the r m in order to find the first two moments of 
/ or S p . The density of time-delays is p3| 



P(r) 



N 

E 

rn— 1 



(S(Trn T)) 



N 



V( T + - T )i T 



t± = 27r7i(3± V&)/NA 

The pair correlation function Ki{t\,t-i) is a universal 
function of the arguments T\ and ti and the "spectrum 
edges" r_ and r + With the help of the pair cor- 

relation function [45| we find that, up to corrections of 
order 1/A 4 , 




A 
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